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We report a Grand Canonical Monte Carlo simulation study of the liquid-vapor interface of a 
model fluid exhibiting polydispersity in terms of the particle size a. The bulk density distribution, 
p^{ff), of the system is controlled by the imposed chemical potential distribution ^{a). We choose 
the latter such that p^{a) assumes a Schulz form with associated degree of polydispersity ~ 14%. 
By introducing a smooth attractive wall, a planar liquid-vapor interface is formed for bulk state 
points within the region of liquid-vapor coexistence. Owing to fractionation, the pure liquid phase is 
enriched in large particles, with respect to the coexisting vapor. We investigate how the spatial non- 
uniformity of the density near the liquid- vapor interface affects the evolution of the local distribution 
of particle sizes between the limiting pure phase forms. We find (as previously predicted by density 
functional theory, Bellier-Castella et al, Phys. Rev. E65, 021503 (2002)) a segregation of smaller 
particles to the interface. The magnitude of this effect is quantified for various a via measurements of 
the relative adsorption. Additionally, we consider the utility of various estimators for the interfacial 
width and highlight the difficulties of isolating the intrinsic contribution of polydispersity to this 
width. 



I. INTRODUCTION AND BACKGROUND 

Complex fluids in which the particles are similar in 
character but not strictly identical, are termed polydis- 
perse. Examples of such arise throughout soft matter 
science, notably in colloidal dispersions, polymer solu- 
tions and liquid-crystals. Typically the polydispersity of 
these systems is manifest as variation in some physical 
attribute such as size, shape or charge, which one de- 
notes by a continuous parameter cr. The form of the 
polydispersity is then quantifiable in terms of a density 
distribution p{a) measuring the number density of par- 
ticles of each a. Accordingly one can regard the system 
as a mixture of an infinite number of particle "species" 
each labelled by the value of a [1] . 

As has long been appreciated, polydispersity can 
deeply influence the thermodynamical and processing 
properties of complex fluids [2-5], making a clear elu- 
cidation of its detailed role a matter of both fundamen- 
tal and practical importance. The majority of recent 
effort in this regard has focussed on clarifying the bulk 
phase behaviour of polydisperse systems (see [6] for a re- 
cent review), which is known to be considerable richer 
in both variety and character than that of correspond- 
ing monodisperse systems. The source of this richness 
can be traced to fractionation effects: at coexistence a 
polydisperse fluid described by some initial "parent" dis- 
tribution, p^{cr), may split into two or more 'daughter' 
phases p'^°'\(t), a = 1, 2, each of which differs in com- 
position from the parent. The sole constraint is that the 
volumetric average of the daughter distributions equals 
the parent distribution. 

The effect of fractionation on phase diagrams can be 
dramatic. For instance, the familiar liquid-vapor bin- 
odal in the density-temperature plane of a monodisperse 
fluid splits into cloud and shadow curves [6], as shown 
schematically in fig. 1(a). These mark, respectively, the 
density of the onset of phase separation and the density of 



the incipient (shadow) phase. The critical point occurs 
neither at the extremum of the cloud nor the shadow 
curve but at their intersection. One implication of this is 
that even at the critical temperature, liquid vapor coexis- 
tence can occur provided the overall parent density is less 
than its critical value. Additional insight into fractiona- 
tion effects can be gleaned from the pressure-temperature 
plane of the phase diagram. In a monodisperse system, 
coexistence occurs along a line in this plane which ter- 
minates at the critical point. However, the introduction 
of polydispersity broadens this line into a region having 
a 'banana'-like shape [7, 8]. The critical point generally 
lies neither at a point of maximum temperature nor max- 
imum pressure on the perimeter of this region (fig. 1(b)). 
Traversing the coexistence region from one pure phase to 
the other (eg. along an isobar or isotherm) corresponds 
to smoothly varying the relative volumes of the system 
occupied by the two coexisting phases, with concomitant 
smooth variation in the forms of the daughter distribu- 
tions. 

In view of the richness of the bulk phase behaviour, it 
is natural to enquire how polydispersity influences fluid 
interfacial properties. Previous work on this subject has 
largely been limited to the study of size-disperse hard 
spheres either at a single hard wall [9-12] or confined to 
a planar slit [13, 14]. These investigations focussed prin- 
cipally on the effect of the wall on the local fluid struc- 
ture and distribution of particle sizes. However, since 
no explicit interparticle attraction was considered, issues 
of wetting and fluid-fluid coexistence did not arise. The 
sole previous study of a fluid-fluid interface in a polydis- 
perse system (of which we are aware) is the density func- 
tional theory study of Bellier-Castella et al [15]. Building 
on previous studies of a homogeneous polydisperse van 
der Waals fluid [8], these authors obtained the density 
proflles for various species across a planar liquid-vapor 
interface. Their calculations indicated a preferential ad- 
sorption of small particles at the interface which they 
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Conventionally the state of such a fluid in the bulk is 
described by a parent density distribution p^{(j), which 
can be written [1]: 



p\a)=n^.f{a) 



(2.2) 



Here no = N/V is the overall particle number density, 
while /((t) is a normalized shape function whose average 
value (T serves to set the scale for all lengths. Since the 
form of /(a) is fixed, the bulk phase diagram is spanned 
by no and the temperature T (cf. fig 1(a)). 

A commonly used measure of the scale of variation in 
the particle diameters is provided by the dimensionless 
degree of polydispersity, defined as the standard devia- 
tion of the parent distribution, normalized by its mean: 



(2.3) 



FIG. 1: Schematic phase diagram as described in the text, 
(a) The no-T plane, (b) The p-T plane. 



In the present work, we have elected to study a parent of 
the Schulz form: 



reported to be broadened with respect to the monodis- 
perse limit. 

In the present work, we have attempted to extend the 
understanding of interfacial behaviour in polydispersc 
fluids by performing detailed Monte Carlo simulations 
of a system of spherical size-disperse particles interact- 
ing via a Lennard- Jones potential. The simulations are 
performed within the grand canonical ensemble and em- 
ploy a chemical potential distribution /i(cr), the form of 
which is chosen such as to yield a bulk density distribu- 
tion /0°((t) having a flxed Schulz form, with associated 
degree of polydispersity w 14%. Such "fixed polydisper- 
sity" corresponds to the experimental situation in, for ex- 
ample, colloidal dispersions or polymer solutions, where 
the form of the parent distribution is prescribed by the 
synthesis process of the particles, and only its scale can 
change depending on the quantity of solvent present. 



II. METHOD 

A. Model and observables 

The model we consider comprises a system of particles 
interacting via an interparticle potential of the Lennard- 
Jones (LJ) form: 
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Here the parameter Z controls the width of the distribu- 
tion and thence the value of 6. We have considered the 
case Z = 50, corresponding to 5 = {Z + 1)~^/^ ~ 0.14. 
The resulting form of /(o") is shown in fig. 2. Note that in 
contrast to, for example, a Gaussian, the Schulz distribu- 
tion vanishes smoothly (has a natural cutoff) as tr ^ 0. 
For the purposes of the MC simulations described below, 
however, one does require an upper cutoff in a, beyond 
which /(ct) is truncated. We set this to be ctc = 1.6. 




FIG. 2: The Schulz parent size distribution f{a) studied in 
this work. 

In general, a polydisperse fluid described by eqs. 2.1 
will exhibit liquid-vapor coexistence within the region of 
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the phase diagram enclosed by the cloud curve (cf. fig. 1). 
For such coexistence states, a liquid-vapor interface can 
be formed by introducing a sufficiently attractive wall in 
the plane 2; = 0. In this work, the particle-wall interac- 
tions are assigned the form 
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, (2.5) 



where iTiu; = {<Ti+d-)/2. The form of this potential derives 
from regarding the wall as a mono-layer of monodisperse 
particles having diameter a and interacting with fluid 
particles via eq. 2.1. 

In order to quantify the resulting interfacial properties, 
we consider the ensemble averaged local density distribu- 
tion p{a, z) at a perpendicular distance z from the wall: 



^0 



p{(j, r) dxdy . 



(2.6) 



We shall also find it useful to define a number of one- 
dimensional profiles which derive from p{(j, z). These are 
the overall density profile 

p{z) = j da p{a,z); (2.7) 
the volume fraction profile 



V{z) = I da p{a, z)a^ ; 



(2.8) 



the local concentration profile: 



(i}(a, z) = 



Pi<^,z) 



(2.9) 



and finally, (in analogy with eq. 2.2), the local normalized 
size distribution: 



Pi<^\z) 
J dap{a\z) 



(2.10) 



The utility of all these quantities will become apparent 
in sec. III. 



B. Simulation methodology 

The grand canonical Monte Carlo algorithm utilized 
to study our model deploys four types of operation: par- 
ticle displacements, deletions, insertions, and resizing. 
The particle diameter a is treated as a continuous vari- 
able throughout. Observables, on the other hand (princi- 
pally the form oi p{a,z)) are accumulated as histograms. 



formed by discretising into bins the permitted ranges of 
a < ac (see sec. II A) and z (i.e. < z < L^). The bin 
widths used were 5a = 0.02 and Sz = 0.02. Further de- 
tails of the implementation can be foimd elsewhere [16]. 

The algorithm requires as input a chemical potential 
distribution fj,{a). Given a nominated shape function 
f{a), specification of the overall density no and the tem- 
perature T serves to fix the bulk form of p,{a\no, T). The 
assumption that the confined fluid exists in equilibrium 
with a bulk reservoir implies that both the bulk and con- 
fined systems have equal n{a). Thus in order to perform 
simulations of the fluid at a wall, we first require the 
form of p,{a) at specified points inside the bulk coexis- 
tence region. For the present work, we have focused on 
the critical isotherm Tc = fc^Tc/e = 1.384, the critical 
temperature having been determined in a separate study 
[17]. Recall (sec. I) that Hquid- vapor coexistence can 
occur in a polydisperse fiuid even at the critical temper- 
ature, provided that the overall density no lies between 
the densities of the cloud point and critical point. At Tc, 
the cloud point occurs at no = 0.044(1) and the critical 
point at no = 0.326(2). 

Using fully periodic simulations, we have determined 
the form of p{a) for a selection of values of no along 
the critical isotherm between the cloud point and the 
critical point. For this purpose, the combined techniques 
of non-equilibrium potential refinement and histogram 
reweighting were utilized. A description of the procedure 
have been previously presented in refs. [18-20] and we 
again refer the interested reader to these papers for full 
details. 

The resulting forms of ii{a\no,Tc) for the respective 
values of no were then used to study the effects of intro- 
ducing two oppositely facing walls at 2; = and z = 
(the system remaining periodic in the x and y directions). 
The system size was set at = Ly = 15a, = 40ct. 
Interactions between the particles and the wall at z = 
were assigned the form eq. 2.5, while a simple hard (im- 
penetrable) wall condition was applied bX, z = L^. This 
arrangement ensures that only the wall aX z = Q can 
become wet. 

For values of no near the cloud point, no liquid-like 

layer was formed at the attractive wall. However, for 
larger no a moderately thick layer was observed. The 
histogram of the interfacial profile p{a, z) presented be- 
low correspond to the choice of parameters no = 0.089 
f = fc = 1.384, with a waU strength e^^/ksT = 6.0. The 
resulting liquid-like layer was found to have a thickness 
of approximately 14ct. This is sufficiently large that the 
liquid- vapor interface can be considered to be essentially 
decoupled from cither wall, a view that was confirmed by 
comparing the liquid and vapor properties on either side 
of the interface with those obtained in the corresponding 
bulk simulations. 

In the course of the simulations, the form of p{a, z) was 
observed to exhibit large, slow fluctuations in which the 
layer thickness varied in the range \2a to 18(t. These fluc- 
tuations complicate the task of accumulating high statis- 
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FIG. 3; Illustration of the procedure for locating the inter- 
face center, as described in the text. Shown is the running 
integral p{a. z')dz' for three example instantaneous density 
profiles that span the typical range of fluctuations. Dashed 
lines denote linear fits to the running integral in the pure liq- 
uid and vapor regions away from the interface. Dotted arrows 
show the position of the interface center determined from the 
intersection of the fits for each profile. 



tics for the intrinsic profile shape because; an ensemble 
average over independent configurations will be consid- 
erably smeared out in the z direction. To circumvent this 
problem we have accumulated a centered profile p{(7, z*), 
whereby the instantaneous profile p{a, z) is first shifted 
with respect to some nominal origin (center) before be- 
ing accumulated in the histogram. The center of p{a, z) 
was itself determined by measuring the running integral 
p{a, z')dz' . This riuantity manifests a smoothed dis- 
continuity at the interface center, the location of which 
can be estimated from the intersection of respective lin- 
ear fits in the pure vapor and liquid regions well away 
from the interface, as illustrated in fig. 3. We note that 
our centering procedure differs to that adopted in some 
other studies of liquid- vapor interfaces (see eg [21, 22]), 
where the interface origin is defined via the instantaneous 
Gibbs dividing surface. Nevertheless, we do not expect 
differences in the precise definition to have repercussions 
for the qualitative features of results. Indeed, as a check 
of our procedure, we have investigated the effect of dis- 
abling particle transfer MC moves (leaving operational 
only particle displacement and resizing moves) . Doing so 
greatly suppresses the fluctuations in the interface thick- 
ness and yields an average profile, the shape of which 
agrees to within error with that obtained from the cen- 
tering technique. 



III. RESULTS 

Before describing our findings for the liquid-vapor in- 
terfacial properties, it is instructive to quantify the na- 
ture of the fractionation in the pure phases away from the 
interface. This is done in fig. 4, which displays the sin- 



gle liquid and vapor (daughter) phase distributions pl((t) 
and /9v((t), together with the Schulz parent distribution 
from which they derive, all at the nominated bulk den- 
sity of no = 0.089. Clearly there is a pronounced segre- 
gation of larger particles to the liquid phase. Indeed, for 
a — ac the density in this phase exceeds that of the vapor 
by a factor of two orders of magnitude-far greater than 
the ratio of the overall number densities in the coexist- 
ing phases (pl/pv ~ 5). One notes further that, as the 
particle size decreases, the density distributions in both 
phases tend to the parent form. The origin of this feature 
is traceable to the fact that the smallest particles inter- 
act only very weakly (cf. eq. 2.1), and consequently their 
density is principally controlled by the imposed chemical 
potential, rather than the number density. 
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FIG. 4: Density distributions p(a) in the coexisting bulk 
phases at no = 0.088, T = 1.384. Also shown is the par- 
ent form p°(cr). Statistical errors are comparable with the 
line widths. 

Turning now to the interfacial properties, fig. 5 
presents the measured form of the centered distribution 
p{a,z*), accumulated as described in sec. II B. The cor- 
responding forms of p{z*) and 77(2;*) (cf. eqs. 2.7 and 
2.8) are given in fig. 6(a). Also shown (fig. 6(b)) is the 
variation of the average particle size a{z*) across the in- 
terface. 

We have attempted to fit p{z*), r)(z*) and a{z) using 
a function of the standard tanh form [23, 24] : 



y{z*\a) = ya -I- yd tanh 



2(4 - 

X 



(3.1) 



where ya = (vl + 2/y)/2 and ya = {yL - yy)/2, with 
yL and yv the appropriate limiting pure phase quan- 
tity. Here Zq denotes the location of the interfacial mid- 
point, while A is a measure of the interfacial half-width. 
It should be noted that the fitting form eq. 3.1 derives 
from mean field theories of coexistence between symmet- 
rical phases [23, 24], and is not necessarily expected to 
hold for real asymmetric fluids which lack particle-hole 
symmetry [25-27]. Nevertheless cq. 3.1 has been shown 
to provide a good description of the liquid-gas interface 
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FIG. 5: Esimtatcs of the profile p[(j, z*) near the liquid-gas 
interface, determined by the procedure described in the text. 

density profile of the monodisperse L J fluid [28] , and in- 
deed we do appear to obtain good fits for all profiles, as 
evidenced in fig. 6. The associated mid-points and inter- 
facial widths are Zq = 0.661, w = 2\ = 9.3(2) for the p{z) 
profile, = 0.142, w = 9.6(2) for r](z), and zj^ = 1.89, 
w = 10.4(3) for a{z). Discrepancies in these estimates 
are presumably traceable to the varying sensitivity of the 
respective observable to the interface proximity. For ex- 
ample, because p{z*) contains no direct information on 
the variation of the concentration across the interface, 
fits to its form possibly constitute a poorer estimator for 
the interfacial width than &(z*). We shall return to this 
point below and in sec. IV. 

It is of interest to compare the interfacial properties 
of our polydisperse fluid with those of the corresponding 
monodisperse system. For reasons detailed in sec. IV, we 
have chosen to perform this comparison not at a given 
temperature and overall bulk number density, but under 
two separate conditions relating to the volume fractions 
of the coexisting pure phases. Firstly we consider, for the 
monodisperse system, the temperature at which the value 
of the relative fluctuation 2{riL — rjv)/{i]L + i]v) equals 
that pertaining to the polydisperse system at the state 
point under consideration (T = 1.384, no = 0.089). This 
occurs for T = 1.080 [29]. The associated monodisperse 
density profile is included in fig. 6(a), together with (in 
the inset) a comparison of the volume fraction profiles of 
both systems. Fitting the monodisperse volume fraction 
profile via the tanh form (eq. 3.1), yields a width estimate 
w = 8.8(2), ie. some 10% smaller than found for the poly- 
disperse profile w = 9.6(2). An alternative scenario com- 
pares the poly- and mono-disperse systems under condi- 
tions of equal volume fraction difference rji^ — rjy (which 
occurs for the monodisperse system at T = 1.048). Un- 
der such conditions the difference in the width estimates 
falls to 5%. We discuss these findings further in sec. IV. 

Given knowledge of the distribution p{a,z*), it is 




FIG. 6: (a) Estimates of the overall number density profile 
p{z*) and volume fraction profile ri{z*). Lines are fits of the 
form eq. 3.1. Also shown is the density profile in the monodis- 
perse limit at r = 1.080, for which the relative volume frac- 
tion fluctuations match those of the polydisperse case (see 
text). The inset shows a comparison of the volume fraction 
profiles under these conditions. Statistical errors axe smaller 
than the symbol sizes, (b) The variation across the interface 
of the average local particle size, (j{z*). The line is a fit of 
the form eq. 3.1. 



straightforward to extract the profiles for given individ- 
ual species i.e. the forms of p{z*\a-). A representative 
selection of these is presented in figs. 7, from which one 
observes that for smaller values of a, a pronomiccd "seg- 
regation" peak occurs in the density near the interface. It 
is interesting to note that for the smallest particle size for 
which reasonable statistics could be obtained (cr = 0.6), 
the densities on either side of the interface barely dif- 
fer, in accord with the coincidence (noted above) of the 
daughter distributions for small a. Notwithstanding this, 
a clear 10% enhancement of the density of this species 
occurs at the interface. That the profile p{z*\a = 0.6) 
couples to changes in the overall density p{z*) but not 
its absolute value provides an indication that the segre- 
gation phenomenon is associated with the surface tension 
of the interface. This point is discussed further in sec. IV. 

The absolute height of the segregation peak is greatest 
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FIG. 7: (a) Estimates of intcrfacial density profiles p{z*\a) 
for various a in the range 0.9 < a < 1.5. The sohd hne is the 
overall density profile p{z). Statistical errors are smaller than 
the symbol sizes, (b)-(e) show the corresponding profiles for 
a selection of cr in the range 0.6 < a < 0.9. The statistical 
quality of the data becomes poor for a < 0.6 due to the low 
parent density (cf. fig. 4). Statistical errors are comparable 
with the spread of data points. 



bility of the segregation peak, which would otherwise be 
overwhelmed by the large relative variation in the con- 
centration profile across the interface, as indeed occurs 
for smaller a. This situation should be contrasted with 
that for p{z*\(7), where the ratio of densities in the liquid 
and vapor become more closely matched as a decreases 
(cf. fig. 7), thereby enhancing the visibility of the segre- 
gation peak. Additionally, it appears likely that any such 
peak is more distinguishable against a monotonically de- 
creasing profile as is the case for the density p{z*\a), than 
for monotonically increasing ones as is the case for the 
concentration profiles ^{z*\a) for a < 1.1 
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for fT « 0.9 — 1.0 and is apparent in all density pro- 
files for a < 1.0. However, the associated segregation 
effect actually extends to at least a = 1.1, as is evi- 
dent from an examination of plots of the local concen- 
tration (j>{a,z*) shown in fig. 8. Specifically, the profile 
(f>{z*\a = 1.1) shows a clear peak at the interface. The 
fact that a peak is visible in this profile, but neither in 
p{z*\a = 1.1) nor the concentration profiles for smaller 
a values, can be explained as follows. For a < 1.1, the 
concentration profiles are monotonically increasing with 
z* and the ratio of concentration in the vapor to that 
in the liquid increases rapidly with decreasing cr. For 
a > 1.1, however, the profiles are monotonically decreas- 
ing and the ratio of concentration in the vapor to that in 
the liquid decreases rapidly with increasing a. Accord- 
ingly, for (7 w 1.1, the concentrations in the two phases 
are closely matched. This matching enhances the visi- 



FIG. 8: A selection of profiles of the local concentration 
(l>{a,z*), as described in the text. Statistical errors do not 
exceed the line widths. 

A more quantitative approach to the segregation effect 
is provided by a quantity known as the symmetrized sur- 
face segregation, first introduced in the context of binary 
liquid mixtures [30, 31]. This compares the intcrfacial 
density profiles of a nominated species a with that of 
some reference species a: 



where the symmetrized concentrations are defined as 
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Here PL(cr) and py(cT) denote the density of a given 
species in the bulk Uquid and vapor phases respectively. 

A(^*|(j, (t) is, by construction, zero in the two pure 
phases, while at any point in the interfacial region it is 
either negative or positive depending on the sign of the 
relative segregation of the two species a and a. The 
relative adsorption of species a with respect to species a 
can then be calculated as [30] : 



T{a\a) 
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O — < 



(3.5) 



The density and concentration profiles discussed above 
show evidence of surface segregation for ct < 1.1. In 
seeking to quantify surface adsorption, it therefore seems 
reasonable to adopt as the reference species a = 1.2, for 
which no segregation was discernible. The corresponding 
forms of A(2;*|cr, (t) are shown in fig. 9(a), together with 
the relative adsorption r(cr|(7) in fig. 9(b). The latter 
figure clearly shows the adsorption peaking at around 
(T = 0.9 — 1.0. The population of particle whose size 
cr > 1.2 are depleted in comparison. We note that the 
adsorption can be readily calculated via alternative ap- 
proaches, such as that based on the construction of a 
Gibbs dividing surface for the reference species [15, 32]. 
As a check on our procedure, we have also applied this 
approach, finding results that are numerically practically 
indistinguishable to those obtaining from eq. 3.5. 

In a previous study of the liquid- vapor interface of a bi- 
nary fluid mixture [31], the forms of A(^*) was argued to 
provide a more reliable measure of the interfacial width 
than that provided by tanh fits to the overall density 
profile p[z*) alone, which tend to underestimate the true 
width. Adopting the criterion [31] that the width is given 
by the range of z* for which |A(z*)| exceeds 5% of its 
maximum value, knowledge of A(z*|a, a) provides esti- 
mates of the interfacial width w(cr) for each u. This 
quantity exhibits a spread of values, the largest corre- 
sponding to the largest species, w{a = 1.6) = 10.4(3), 
a values which indeed clearly exceeds that of w = 9.3(2) 
arising from the tanh fit to the density profile p{z*). It is, 
however, somewhat closer to that obtained from the fit to 
ri(z*) (i.e. u] = 9.6(2)), and agrees well with the estimate 
w = 10.4(2) obtained for the fit to &{z*). We postpone 
further discussion of this comparison until sec. IV. 

Having examined the density and concentration pro- 
files and associated adsorption, we turn finally to consider 
the variation of the normalized size distribution /(u) 
across the interface, eq. 2.10. The forms of f{a\z*) for 
selected values of z* spanning the interface are shown in 
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FIG. 9: (a) Estimates of the symmetrized surface segregation 
A(2*|a, (t) (see text) for a choice of the reference species a = 
1.2. Statistical errors are comparable with the symbol size, 
(b). The cr-dependent adsorption relative to a = 1.2. 



fig. 10(a). From the figure, one observes that on travers- 
ing from the pure liquid to the pure vapor, the size dis- 
tribution steadily narrows, while its mean shifts to lower 
a (cf. figs. 4, 6(b)). The scale of the evolution of the 
fractionation across the interface becomes more apparent 
when referenced with respect to one of the pure phases 
(we have chosen the vapor) and plotted on a log scale, as 
shown in fig. 10(b). 



IV. DISCUSSION AND CONCLUSIONS 

In summary, we have performed a detailed simulation 
study of the liquid- vapor interface of a polydisperse fluid. 
The main feature of our results is the finding of a prefer- 
ential adsorption of smaller particle at the interface and 
a corresponding depletion of larger ones. These findings 
are in broad accord with those of a density functional 
theory for a polydisperse van-der-Waals fluid [15]. We 
note however that the segregation phenomenon is not 
peculiar to the case of polydisperse mixtures. Analogous 
effects are well known to occur in the context of binary 
fluid mixtures, where one component is often preferen- 
tially absorbed at the liquid-gas boundary (see eg. refs. 
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FIG. 10: (a) Variation of the normalized local size distri- 
bution f{a\z*) across the interface, as described in the text. 
Lines through symbols are merely guides to the eye. (b) The 
same data referenced with respect to the vapor phase distri- 
bution, and expressed on a log scale. 

[22, 30, 33-36]). The segregation occurs because coating 
the interface in the more volatile (less strongly interact- 
ing) of the two species reduces the surface tension, as 
can be shown from the Gibbs adsorption equation [32]. 
In view of this, one might expect that for the polydis- 
perse system the smallest species should be maximally 
adsorbed. By contrast, we find that the adsorption is ac- 
tually greatest for an intermediate size (cf. fig. 9). This 
finding can be understood from the fact that for our inter- 
particle potential (eq. 2.1) the smallest particles interact 
only weakly with one another (and indeed become ideal 
in the limit a ^ 0). Coating the interface in such very 
small particles, cannot therefore screen the larger ones 
and reduce the surface tension. 

We have examined various quantities which provide a 
potential measure of the liquid-vapor interfacial width. 
Good fits of the tanh form (eq. 3.1) were achieved for 
the profiles of the overall number density, volume frac- 
tion and local average particle size. Additionally, a fur- 
ther measure, the 95% width of the symmetrized surface 
segregation profile, A(z*|(t, (t), was examined [31]. This 
latter quantity contains explicit information on the den- 
sity profiles of the individual species p{z*\a), which, it 
seems, is important in obtaining reliable estimates of the 
interfacial width. By contrast, the overall profiles p{z*) 



and ri(z*), integrate out some, or all. of the a depen- 
dence of p{a, z*). They thus suppress information about 
species whose individual profile width is large, but which 
contribute relatively little weight to the overall density. 
This may result in an underestimate of the true length 
scale over which the system properties deviate from their 
bulk values, and indeed the discrepancies we find in the 
interfacial width for the various profiles (sec. Ill) provide 
some evidence for this. We note further that if one sim- 
ply examines the overall density profile, one may miss 
important features such as the segregation peak, which 
is apparent in p{z*\a) for small cr, but not in p{z*). 

Turning now to the effects of polydispersity on the 
width of the liquid- vapor interface, it is stated in ref.[15] 
that polydispersity broadens the interface compared to 
the monodisperse limit. However, our results highlight 
the need for care when attempting to isolate the in- 
trinsic effects of polydispersity on the interfacial width 
from those arising indirectly as a result of polydispersity- 
indiiccd alterations to the bulk phase behaviour. In- 
deed, the latter may well constitute the dominant fac- 
tor. For example, introducing polydispersity of the form 
described by eq. 2.1 (which is similar to that considered 
in [15]) tends to raise the critical temperature sharply 
with respect to the monodisperse limit [8, 17]. Thus on 
increasing the degree of polydispersity at a given tem- 
perature, one can expect the interface to become sharper 
since the system moves deeper into the two phase region. 
Our results provide an extreme example of this. Specifi- 
cally the temperature at which we studied phase coexis- 
tence in the polydispersc system, T = 1.384 exceeds the 
corresponding critical temperature in the monodisperse 
limit, namely = 1.1876 [29]. Thus at T = 1.384 no 
interface would occur at all in the monodisperse limit. 

Another issue which complicates comparison of inter- 
facial widths for polydisperse and monodisperse system 
is the choice of parent density uq. Owing to fraction- 
ation and the resulting smearing of the phase diagram 
(cf. figs. 1), the value of no represents a crucial fac- 
tor in determining the properties of the coexisting bulk 
phases and hence the interface between them. Again the 
present work provides a case in point: although we stud- 
ied the system at the critical temperature, true phase 
coexistence was nevertheless observable because hq was 
set less than its critical value. Indeed one can obtain a 
sharper or a broader interface simply by changing no at 
constant T, a situation which contrasts with that in the 
monodisperse limit. 

Thus we believe that when trying to isolate the in- 
trinsic effect of polydispersity on interfacial widths, one 
should endeavor to perform the comparison at points in 
the phase diagram that are in some sense equivalent in 
terms of their relative deviation from criticality. Match- 
ing the bulk volume fraction fluctuations (in preference 
to the densities which do not provide equivalent infor- 
mation) seems one reasonable way to do this. Doing so 
(cf. inset of fig. 6(a)) shows that polydispersity may only 
broaden the interface marginally (5% — 10%). We cau- 
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tion, however, that our study is far from comprehensive 
in this regard, and it is possible that interfacial segre- 
gation and its effect on the surface tension and interfa- 
cial width might become more significant at state points 
further removed from criticality [15]. This would be an 
interesting avenue for further computational study. 
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